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Abstract  Multiple  oximes  have  been  synthesized  and 
evaluated  for  use  as  countermeasures  against  chemical 
warfare  nerve  agents.  The  current  U.S.  military  and  civilian 
oxime  countermeasure,  2-[(hydroxyimino)methyl]-l-meth- 
ylpyridin-l-ium  chloride  (2-PAM),  is  under  consideration 
for  replacement  with  a  more  effective  acetylcholinesterase 
reactivator,  1,1’  -methylenebis  { 4-hydroxyiminomethyl  }py- 
ridinium  dimethanesulfonate  (MMB-4).  Kinetic  data  in  the 
scientific  literature  for  MMB-4  are  limited;  therefore,  a 
physiologically  based  pharmacokinetic  (PBPK)  model  was 
developed  for  a  structurally  related  oxime,  l,l’-trimethyl- 
enebis{4-hydroximinomethyl}pyridinium  dibromide.  Based 
on  a  previous  model  structure  for  the  organophosphate  dii- 
sopropylfluorophosphate,  the  model  includes  key  sites  of 
acetylcholinesterase  inhibition  (brain  and  diaphragm),  as 
well  as  fat,  kidney,  liver,  rapidly  perfused  tissues  and  slowly 
perfused  tissues.  All  tissue  compartments  are  diffusion  lim¬ 
ited.  Model  parameters  were  collected  from  the  literature, 
predicted  using  quantitative  structure-property  relationships 
or,  when  necessary,  fit  to  available  pharmacokinetic  data 
from  the  literature.  The  model  was  parameterized  using  rat 
plasma,  tissue  and  urine  time  course  data  from  intramuscular 
administration,  as  well  as  human  blood  and  urine  data  from 
intravenous  and  intramuscular  administration;  sensitivity 
analyses  were  performed.  The  PBPK  model  successfully 
simulates  rat  and  human  data  sets  and  has  been  evaluated  by 


T.  R.  Sterner  (El)  •  C.  D.  Ruark  •  T.  R.  Covington  • 

J.  M.  Gearhart 

Henry  M.  Jackson  Foundation  for  the  Advancement 
of  Military  Medicine  (HJF),  2729  R  Street,  Bldg  837, 
Wright-Patterson  AFB,  OH  45433-5707,  USA 
e-mail:  teri.stemer.ctr@wpafb.af.mil 

K.  O.  Yu 

AFRL  711  HPW/RHDJ,  Wright-Patterson  AFB,  OH,  USA 


predicting  intravenous  mouse  and  intramuscular  human  data 
not  used  in  the  development  of  the  model.  Monte  Carlo 
analyses  were  performed  to  quantify  human  population 
kinetic  variability  in  the  human  evaluation  data  set.  The 
model  identifies  potential  pharmacokinetic  differences 
between  rodents  and  humans,  indicated  by  differences  in 
model  parameters  between  species.  The  PBPK  model  can  be 
used  to  optimize  the  dosing  regimen  to  improve  oxime 
therapeutic  efficacy  in  a  human  population. 

Keywords  Oxime  •  PBPK  model  •  Rat  •  Human  ■ 
Pharmacokinetic  •  Organophosphate 

Introduction 

Traditional  oximes  are  pyridinium  mono-  or  bis-quaternary 
compounds  characterized  by  their  quaternary  ammonium 
structures  and  low  lipophilicity  (negative  log  P  values) 
(Voicu  et  al.  2010).  Medically,  oximes  are  administered  to 
counteract  organophosphate  (OP)  poisoning.  OPs  form 
serine-conjugated  phosphonates  or  phosphorates  with  ace¬ 
tylcholinesterase  (AChE).  Oxime  molecules  bind  to 
OP-AChE  complexes  and  regenerate  AChE  activity 
through  nucleophilic  attack  of  the  oximate  on  the  phos¬ 
phorus  atom,  leading  to  a  pentacoordinate  phospho- 
ryloxime-enzyme  transition  state.  The  unstable  geometry 
of  this  complex  liberates  the  enzyme  from  the  phospho- 
ryloxime  (Eyer  2003).  The  resulting  reactivation  of  AChE 
enzymes  restores  cholinergic  function  disrupted  by  ace¬ 
tylcholine  buildup  in  the  post-synaptic  clefts  of  neuro¬ 
muscular  junctions  (Saxena  et  al.  2008). 

While  used  in  clinical  practice  as  part  of  the  treatment 
for  OP  pesticide  exposure  (Eddleston  et  al.  2008),  effective 
oxime  treatment  is  also  a  major  concern  for  military 
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members  in  the  event  of  a  chemical  warfare  nerve  agent 
(CWNA)  attack.  2-[(Hydroxyimino)methyl]-l-methylpyri- 
din-l-ium  chloride  (2-PAM)  has  been  fielded  with  U.S. 
troops  since  the  1970s  (Saxena  et  al.  2008).  The  MARK  I 
field  kit,  used  by  U.S.  military  and  civilian  agencies,  is  a 
double-barreled  autoinjector  containing  a  2-mg  atropine 
sulfate  dose  alongside  a  600-mg  dose  of  2-PAM.  Current 
military  protocol  is  to  use  three  MARK  I  kits  as  a  loading 
dose  in  a  severely  affected  victim,  without  exceeding 
1,800  mg  2-PAM  per  hour.  Also  in  the  kit  is  a  separate 
autoinjector  containing  10  mg  diazepam  as  an  antidote  for 
convulsions.  Intramuscular  autoinjectors  allow  for  rapid 
self-  or  buddy-administration  into  a  large  muscle  such  as 
the  thigh  (Newmark  2009).  Traditional  oximes  must  be 
administered  intravenously  or  intramuscularly  because 
their  oral  uptake  is  poor,  as  indicated  by  the  negative  log  P 
values.  Oral  bioavailability  of  several  oximes,  including 
2-PAM,  is  only  about  6  %  of  the  intramuscular  bioavail¬ 
ability  (Voicu  et  al.  2010). 

2-PAM  is  under  consideration  for  replacement  with  a 
more  effective  AChE  reactivator  since  the  efficacy  of 
2-PAM  is  extremely  limited  against  CWNAs  such  as  soman 
(GD)  and  cyclosarin  (GF)  (Lundy  et  al.  2011;  Saxena  et  al. 
2008).  A  strong  candidate  for  replacement  is  the  symmetric 
bispyridinium  oxime  1 , 1  ’-methylenebis  { 4-hydroxyi  mi  no- 
methyl  jpyridinium  dimethanesulfonate  (methoxime  or 
MMB-4)  (Lundy  et  al.  201 1).  MMB-4  is  effective  against  a 
broad  range  of  OP  CWNAs,  including  GD  and  GF,  as  well  as 
other  well-known  agents  such  as  tabun  (GA),  sarin  (GB)  and 
Russian  VX  (VR)  (Koplovitz  et  al.  2007;  Saxena  et  al.  2008). 
MMB-4  outperformed  eight  other  oximes,  including 
2-PAM,  in  reactivating  GB,  GF,  VX  and  VR-inhibited  gui¬ 
nea  pig  AChE  in  blood,  brain  and  three  peripheral  tissues 
(diaphragm,  heart  and  skeletal  muscle)  (Shih  et  al.  2009). 
MMB-4,  having  a  higher  LD50  in  both  rats  and  mice  (Eyer 
and  Worek  2007),  is  also  less  toxic  than  2-PAM  in  these 
rodents.  While  MMB-4  is  currently  fielded  in  the  Czech 
Republic,  it  has  not  been  approved  for  use  clinically  in  the 
United  States  (Lundy  et  al.  2011).  For  the  U.S.  military  to 
deploy  an  antidote,  it  must  undergo  full  approval  by  the  U.S. 
Food  and  Drug  Administration  (FDA)  (Newmark  2009). 

Physiologically  based  pharmacokinetic  models  (PBPK) 
are  uniquely  suited  for  extrapolating  kinetics  between  spe¬ 
cies,  routes  of  administration  and  doses  (Loizou  et  al.  2008), 
are  essential  for  integrating  in  vitro  animal  data  to  optimize 
therapeutic  dosing  regimens  in  humans,  and  have  even  been 
used  as  a  tool  to  facilitate  FDA  decisions  (Zhao  et  al.  2011). 
Thus,  in  order  to  better  understand  potential  species-  and 
route-specific  differences,  a  PBPK  model  for  oximes  was 
developed;  however,  development  of  a  PBPK  model  requires 
pharmacokinetic  data.  Kinetic  data  for  MMB-4  are  limited 
at  this  time,  but  sufficient  data  are  available  for  a  closely 
related  oxime,  1,1’  -trimethylenebis  { 4-hydroximinomethyl } 


MMB-4 


TMB-4 

Fig.  1  Chemical  structures  of  MMB-4  and  TMB-4 

pyridinium  dibromide  (trimedoxime  or  TMB-4).  TMB-4, 
while  not  as  pharmacologically  effective  as  MMB-4  against 
select  agents  (GB,  GF,  VX,  VR)  (Shih  et  al.  2009),  is  capable 
of  activating  OP-inhibited  AChE  and  is  actually  more 
effective  against  GA  in  a  human  erythrocyte  in  vitro  study 
(Kovarik  et  al.  2007).  TMB-4  and  MMB-4  are  structurally 
similar  (Fig.  1);  therefore,  in  the  absence  of  sufficient  MMB- 
4  data,  the  oxime  model  was  developed  specifically  for 
TMB-4  as  a  surrogate. 

Developed  herein  is  a  PBPK  model  for  TMB-4,  using  rat 
and  human  datasets  from  the  literature.  The  model  is 
evaluated  against  intravenous  mouse  and  intramuscular 
human  data  not  used  in  the  original  development  of  the 
model.  Monte  Carlo  analyses  are  then  utilized  to  quantify 
human  population  kinetic  variability.  Finally,  sensitivity 
analyses  are  performed.  This  TMB-4  PBPK  model  can  be 
expanded  to  include  pharmacodynamics  (PD),  which 
would  describe  the  interaction  of  OPs,  AChE  and  oximes  at 
the  molecular  level. 


Methods 


MMB-4  and  TMB-4  are  water-soluble  quaternary  ammo¬ 
nium  salts.  TMB-4  is  a  kinetic  surrogate  for  MMB-4  based 
not  only  on  structure  (Fig.  1),  but  also  on  predicted  dis¬ 
tribution  parameters  (Table  1).  Log  P,  pKa  and  log  D 
values  predicted  using  ChemSketch  (Advanced  Chemistry 
Development,  Inc.,  Toronto,  Ontario,  Canada)  indicate 
high  water  solubility  for  both  oximes.  This  property  results 
in  nearly  identical  predicted  tissue/blood  partition  coeffi¬ 
cients  (PCs)  for  the  two  oximes  (Table  1). 

PCs  are  estimated  using  a  modification  of  the  quantita¬ 
tive  structure-property  relationship  (QSPR)  algorithm  by 
Schmitt  (2008)  (Eq.  1): 


K,:pl  = 


Ct_  =  (Q\  (CA 

CPi  \cj  \Cpi) 

( F int Cint  T  FCell(-'cell'\ 

V  ^  j 


FCell\ 

F7 


(/f) 
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where  C,  is  the  chemical  concentration  in  the  tissue,  Cpj  is 
the  chemical  concentration  in  the  plasma,  C„  is  the 
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Table  1  Predicted  distribution  parameters  for  MMB-4  and  TMB-4 


Parameter 

MMB-4 

TMB-4 

Log  Pa 

-1.87 

-3.20 

pKaa 

7.79 

8.03 

Log  Da  b 

-1.22  to 

-1.13 

-2.76  to  -2.68 

Tissue/blood 

partition 

coefficients 

Rat 

Human 

Rat 

Human 

Brain 

0.91 

0.92 

0.91 

0.92 

Fat 

0.27 

0.24 

0.21 

0.19 

Kidney 

0.85 

0.94 

0.84 

0.94 

Liver 

0.69 

0.90 

0.69 

0.90 

Muscle 

0.89 

0.94 

0.89 

0.94 

a  As  predicted  by  ChemSketch 
b  Normal  human  blood  pH  is  7.32  to  7.43  (Wu  2006) 


unbound  chemical  concentration,  Fint  is  the  fraction  of  tissue 
that  is  interstitium,  Cint  is  the  chemical  concentration  in  the 
interstitium,  FcM  is  the  fraction  of  tissue  that  is  cellular,  Cceu 
is  the  chemical  concentration  in  the  cellular  space,//"'  is  the 
fraction  of  unbound  chemical  in  the  interstitium,  '  is  the 
fraction  of  unbound  chemical  in  the  cell  and  //'  is  the  fraction 
of  unbound  chemical  in  the  plasma. 

The  modification  calculates  the  unbound  chemical 
fraction,  /„,  by  accounting  for  the  fraction  of  water  in  each 
respective  tissue,  Fw,  based  on  the  assumption  that  the 
unbound  chemical  is  in  the  water  fraction.  Substituting  into 
the  Schmitt  (2008)  equation  (Eq.  1)  results  in  Eq.  2: 


K,.Pi  = 


FintCir 


F cellC Cell\ 


Cu 


' cell 


j-int  p 

T7  w  _i_  J7  w 

r‘nt  fint  ^  rceU  fcell 

J  u  J  it 


Cu 

Cpi 

ft 

7 pl 


(2) 


where  F™*  is  the  fraction  of  water  in  the  interstitium,  F"‘n 
is  the  fraction  of  water  in  the  cellular  space  and  Frf  is  the 
fraction  of  water  in  the  plasma.  To  convert  Kt:p[  to  a  steady 
state  tissue/blood  partition  coefficient  (Kt:b),  a  red  blood 
cell  (RBC)  compartment  correction  is  made  (Eq.  3): 


Kt,b  = 


Kpp, 


(1  —  Frbc)  +  {K RHO.pl  X  Frbc ) 


(3) 


where  FRBC  is  the  fraction  of  RBCs  in  blood,  as  determined 
from  the  hematocrit  measurements,  and  KRBC:pi  is  the 
steady  state  RBC/plasma  partition  coefficient.  KRBC;p/  may 
be  calculated  using  Eq.  2  where  Fint  is  zero  and  Fcen  is  one. 


PBPK  model  structure 


PBPK  models  simulate  the  absorption,  distribution,  metab¬ 
olism  and  excretion  of  chemicals.  Groups  of  equations  rep¬ 
resent  organs  or  grouped  tissues  that  affect  the 


A 


B 


— 1 

> 

Tissue  Blood 

^ _ 

A 

'  Tissue 

Fig.  2  Schematic  of  the  PBPK  model  structure  for  the  oxime  TMB-4. 
Tissue  compartments  depicted  in  Schematic  a  include  diffusion 
limitation  as  illustrated  by  the  generic  tissue  in  Schematic  b 

pharmacokinetics  of  the  chemical.  The  PBPK  model  pre¬ 
sented  here  (Fig.  2)  is  based  upon  a  previously  published 
PBPK  model  for  an  OP,  diisopropylfluorophosphate  (DFP) 
(Gearhart  et  al.  1990, 1994).  The  model  retains  likeness  to  the 
DFP  model  in  order  to  allow  for  future  prediction  of  phar¬ 
macokinetic/pharmacodynamic  interactions  between  OPs, 
oximes  and  AChE.  The  separate  tissue  compartments  are 
described  by  a  system  of  ordinary  differential  equations  that 
are  solved  simultaneously  using  Gear’s  stiff  algorithm 
implemented  in  acslX  (AEgis  Technologies,  Huntsville,  AL). 

The  TMB-4  model  structure  (Fig.  2)  includes  brain  and 
diaphragm  compartments,  both  of  which  are  OP  target 
tissues  (Gearhart  et  al.  1990,  1994),  and  therefore,  targets 
for  oxime  reactivation.  The  model  also  includes  compart¬ 
ments  for  fat,  kidney  and  liver,  with  the  remaining  organs 
lumped  into  rapidly  and  slowly  perfused  tissue  compart¬ 
ments.  Kidney  and  liver  compartments  represent  not  only 
locations  for  clearance  of  oximes,  but  are  also  tissues  for 
which  oxime  time  course  data  exist.  Fat  was  kept  separate 
from  the  slowly  perfused  tissues  due  to  the  extreme  water 
solubility  of  oximes  which  would  largely  prevent  their 
uptake  into  fat,  but  would  not  inhibit  oximes  from  being 
taken  up  into  more  hydrophilic  slowly  perfused  tissues 
such  as  muscle.  Each  compartment  or  tissue  in  the  model 
(Fig.  2a)  is  coded  with  diffusion-limited  distribution 
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(Fig.  2b)  dependent  not  only  on  the  tissue/blood  PC,  which 
regulates  the  concentration  gradient-driven  diffusion 
across  the  membranes,  but  also  on  the  permeability  area 
(PA)  cross  product  for  the  tissue,  which  slows  the  pene¬ 
tration  of  molecules  that  do  not  readily  diffuse. 

The  model  uses  Michaelis-Menten  kinetics  to  describe 
both  metabolism  in  the  liver  and  excretion  from  the  kidney. 
Minor  amounts  of  TMB-4  metabolism  have  been  shown  by 
de  Miranda  et  al.  (1972)  and  Morgan  et  al.  (1982).  The 
model  also  incorporates  saturable  elimination  kinetics  as 
oximes  are  dependent  on  organic  cation  transporter  (OCT) 
uptake  in  the  proximal  renal  tubules  to  facilitate  excretion 
(Shu  2011). 

The  model  allows  for  TMB-4  dosing  from  intravenous 
and  intramuscular  routes.  Intravenous  injections  are  mod¬ 
eled  as  short  infusions  directly  into  the  venous  blood 
supply.  The  length  of  the  injection  is  set  to  the  study’s 
actual  infusion  time,  if  reported,  or  fit  to  the  data.  A  sub¬ 
compartment  for  intramuscular  dosing  (IM  compartment) 
is  separated  from  the  slowly  perfused  tissue  compartment 
as  needed.  The  model  simulates  a  concentrated  amount  of 
TMB-4  added  quickly  to  a  small  portion  of  muscle.  TMB-4 
then  diffuses  into  the  venous  blood  directly,  rather  than 
coming  to  equilibrium  in  the  entire  slowly  perfused  com¬ 
partment  before  diffusing  into  the  blood.  The  length  of  the 
intramuscular  injection  is  also  set  to  study  specifications, 
when  published,  or  fit  to  the  data.  The  volume  of  the  IM 
compartment  is  calculated  from  the  volume  of  the  intra¬ 
muscular  injection  in  each  study.  This  built-in  flexibility 
allows  the  model  to  adapt  to  different  dosing  protocols.  For 
intravenous  dosing,  the  parameters  for  the  IM  compartment 
are  set  to  zero,  and  the  IM  compartment  is  not  separated 
from  the  slowly  perfused  tissue  compartment. 

TMB-4  is  available  in  solid  forms  as  the  dibromide 
(Milic  et  al.  1996)  and  dichloride  (Calesnick  and  Chris¬ 
tensen  1967;  Vojvodic  1970).  Doses  in  the  TMB-4  studies 
were  reported  as  a  dose  of  a  specific  salt  or  of  the  cation 
only.  Therefore,  all  doses  were  converted  to  mg/kg  body- 
weight  (BW)  of  cation  only,  as  needed,  for  use  in  the 
model.  The  oxime  TMB-4  has  a  molecular  weight  of 
286.33  g/mol  without  its  anion.  For  simplicity,  the  model 
does  not  accommodate  the  different  ionization  states. 

Model  parameterization 

Physiological  parameter  values  (Table  2)  for  each  tissue 
were  found  in  the  published  literature  for  each  species  (rat, 
mouse,  human)  whenever  possible.  Exceptions  and  calcu¬ 
lated  values  are  noted.  When  literature  values  were 
unavailable,  parameters  were  fit  using  the  model  and  the 
model  parameterization  data  sets.  All  simulations  use 
study-specific  BW  values.  Allometric  scaling  (i.e.,  volumes 
scaled  by  BW;  flows  and  metabolic  capacities  scaled  by 


Table  2  Physiological  parameters  used  in  the  TMB-4  model 


Parameter 

Mouse 

Rat 

Human 

Fraction  hematocrit  in 

HEM 

0.45a 

0.46a 

0.44a 

blood 

Flow  rates  (L/h/kg0  75) 
Cardiac  output 

QCC 

16.22b 

14.09b 

12.89 

Alveolar  ventilation 

QPC 

29.091 

19.856b 

12.396b 

Blood  flow  rates  (fraction  of  cardiac  output) 

Brain 

QBrC 

0.033 

0.020 

0.114 

Diaphragm 

QDC 

0.006c 

0.006° 

0.006d 

Fat 

QFC 

0.070e 

0.070 

0.052 

Kidneys 

QKC 

0.091 

0.141 

0.175 

Liver 

QLC 

0.161 

0.183 

0.227 

Rapidly  perfused* 

QRC 

0.59f 

0.47f 

0.66f 

Slowly  perfused* 

QSC 

0.4  lf 

0.53f 

0.34f 

Tissue  volume  (fraction  of  body  weight) 

Arterial  blood 

VAC 

0.0122 

0.0185 

0.0308 

Venous  blood 

VVC 

0.0368 

0.0555 

0.0482 

Brain 

VBrC 

0.017 

0.006 

0.020 

Diaphragm 

VDC 

0.003° 

0.003° 

0.003d 

Fat 

VFC 

0.070 

0.070 

0.214 

Kidneys 

VKC 

0.017 

0.007 

0.004 

Liver 

VLC 

0.055 

0.034 

0.026 

Rapidly  perfused* 

VRC 

0.101f 

0.055f 

0.063f 

Slowly  perfused* 

VSC 

0.793f 

0.821f 

0.844f 

Tissue  blood  volume  (fraction  of  tissue  volume) 

Brain 

VBrBC 

0.03 

0.03 

0.04 

Diaphragm 

VDBC 

0.04s 

0.04s 

0.01s 

Fat 

VFBC 

0.0 15h 

0.015b 

0.020 

Kidneys 

VKBC 

0.24 

0.16 

0.36 

Liver 

VLBC 

0.31 

0.21 

0.11 

Rapidly  perfused 

VRBC 

0.208* 

0.208* 

0.170* 

Slowly  perfused 

VSBC 

0.033* 

0.033* 

0.038* 

IM  compartment  parameters 

Blood  flow  rate  (L/h/g) 

QIMg 

NA-iv 

only 

0.012J 

0.045 

Injection  volume  (mL) 

VolIM 

NA-iv 

only 

0.25k 

2.00* 

Tissue  blood  volume 

VIMBC 

NA-iv 

0.04h 

0.01s 

(fraction) 

only 

All  values  are  from  Brown  et  al.  (1997)  unless  otherwise  noted 
a  Davies  and  Morris  (1993);  b  Arms  and  Travis  (1988);  0  Gearhart 
et  al.  (1990);  d  Gearhart  et  al.  (1994);  °  Rat  value  from  Brown  et  al. 
(1997);  f  Sum  of  applicable  tissue  values  published  in  Brown  et  al. 
(1997);  g  Muscle  value  from  Brown  et  al.  (1997);  h  Guinea  pig  value- 
Bosse  and  Wassermann  (1970);  1  Average  of  applicable  tissue  values 
published  in  Brown  et  al.  (1997); j  Copp  et  al.  (2010);  k  OCV  (201 1); 
1  Vojvodic  (1970) 

*  Total  flow  or  volume  (i.e.,  includes  values  for  brain,  fat,  liver,  etc.); 
the  model  subtracts  the  appropriate  blood  flows  and  tissue  volumes 
from  each  lumped  compartment  following  allometric  scaling 
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BW  to  the  three-quarters  power)  coded  within  the  model 
was  used  to  modify  the  inputs  to  accommodate  interspecies 
extrapolation  (Gearhart  et  al.  1990). 

The  IM  compartment  values  in  Table  2  are  presented  in 
different  units  due  to  the  need  to  calculate  the  size  and  flow 
of  the  compartment  based  on  the  volume  injected  intra¬ 
muscularly  in  each  study.  The  blood  flow  to  the  rat  IM 
compartment  is  set  to  a  rate  of  0.012  L/h/g  tissue  measured 
in  the  thigh  muscle  of  control  group  Wistar  rats  (Copp  et  al. 
2010).  The  fraction  of  cardiac  output  to  the  IM  compart¬ 
ment  is  then  calculated  within  the  code  based  on  the  vol¬ 
ume  of  the  injection.  The  human  muscle  blood  flow  rate  of 
0.045  L/h/g  muscle  tissue  is  calculated  from  the  fraction  of 
cardiac  output  flow  to  muscle,  the  cardiac  output  flow  rate 
and  the  fraction  of  BW  that  is  muscle  (Brown  et  al.  1997). 
The  IM  compartment  volume  is  calculated  in  the  code  to  be 
twice  the  volume  of  the  injection;  of  this  total  volume,  half 
is  assumed  to  be  muscle  tissue  and  half  is  the  injection 
solution.  Although  a  more  physiologically  correct  com¬ 
partment  would  increase  and  decrease  the  volume  with 
injection  and  subsequent  absorption,  this  simplification 
helps  capture  the  effect  of  the  muscle  fibers  being  swelled 
or  separated  to  twice  their  normal  size  while  retaining  their 
original  amount  of  blood  flow. 

Physicochemical  parameters  used  in  the  model  are  sum¬ 
marized  in  Table  3.  PCs  for  TMB-4  were  initially  set  to 
QSPR  estimates  and  were  then  refined  as  needed  by  fitting  to 
the  data.  PCs  were  predicted  for  a  variety  of  tissues  for  both 
rat  and  human  (Table  2);  rat  PCs  are  also  used  for  the  mouse. 
Lumped  compartment/blood  PCs  were  calculated  as  averages 
of  predicted  values  for  rapidly  (heart,  lung)  or  slowly  (mus¬ 
cle,  skin,  bone)  perfused  tissues.  Diffusion  limitation  con¬ 
stants  (permeability  area  cross  products  or  PAs),  metabolism 
and  clearance  parameters  were  fit  to  best  describe  the  rat 
intramuscular  data  set  (Garrigue  et  al.  1991). 

Data  used  in  model  parameterization 

Study  data  for  parameterizing  and  evaluating  the  model  are 
from  the  literature.  A  study  was  utilized  if  data  were  pre¬ 
sented  for  blood  or  plasma  time  course  from  an  intravenous 
or  intramuscular  dose  of  TMB-4  in  any  species.  Only  data 
where  TMB-4  was  administered  alone  were  considered 
since  it  was  assumed  that  co-administration  with  atropine 
would  change  the  kinetics  of  the  oxime  due  to  its  effects  on 
the  physiology  of  the  animal  (e.g.,  cardiac  output). 

Data  from  Garrigue  et  al.  (1991)  were  selected  to 
parameterize  this  PBPK  model  because  the  study  measured 
time  courses  for  plasma,  urine  and  three  tissues.  In  this 
study,  male  Wistar  rats  were  administered  a  single  intra¬ 
muscular  dose  of  14C-labeled  TMB-4  equivalent  to 
21.5  mg/kg  BW.  Rats  weighed  180  to  220  g;  a  midpoint  of 
200  g  BW  was  utilized  for  model  simulations.  Animals 


Table  3  Physicochemical  parameters  used  in  the  TMB-4  model 


Parameter 

Rat/ 

Mouse 

Human 

Metabolism  parameters 

Maximum  enzymatic  metabolic  rate 

VMaxC 

o 

=tfc 

0.6* 

constant  (mg/h/kg0'75) 

Enzyme  affinity  constant  (mg/L) 

Km 

0. 1# 

0.1* 

Urinary  excretion  parameters 

Maximum  saturable  excretion  rate 

VMaxUrC 

80# 

480# 

constant  (mg/h/kg0'75) 

Excretion  affinity  constant  (mg/L) 

KmUr 

900# 

900* 

Dosing  parameters 

Length  of  intramuscular  injection 

TIM 

0.025# 

0.0042“ 

(hours) 

0.0083b 

Length  of  intravenous  injection 

TINF 

0.01c 

0.0625d 

(hours) 

Tissue/blood  partition  coefficients  (unitless f 

Brain 

PBr 

0.1 3# 

0.13* 

Diaphragm 

PD 

0.89 

0.94 

Fat 

PF 

0.21 

0.19 

IM  compartment 

PIM 

7.1# 

2.5* 

Kidney 

PK 

100# 

0.94 

Liver 

PL 

12# 

0.90 

Rapidly  perfused 

PR 

0.73 

0.97 

Slowly  perfused 

PS 

10# 

0.74 

Diffusion  limitation  coefficient  constants  (L/h/kg0'75) 

Brain 

PABrC 

0.002# 

0.002* 

Diaphragm 

PADC 

0.2f 

0.2* 

Fat 

PAFC 

0.38 

0.3* 

IM  compartment 

PAIMC 

12# 

12* 

Kidney 

PAKC 

1.5* 

1.5* 

Liver 

PALC 

0.2* 

0.2* 

Rapidly  perfused 

PARC 

1.5h 

1.5* 

Slowly  perfused 

PASC 

03* 

0.3* 

a  Fit  for  Vojvodic  (1970);  b  Twice  the  value  fit  for  Vojvodic  (1970), 
corresponding  to  two  injections  in  Calesnick  and  Christensen  (1967); 
c  Fit  for  mouse  simulation  (Milic  et  al.  1996);  d  Calesnick  and 
Christensen  (1967)  minimum  study  value;  e  PC  values  predicted 
unless  otherwise  noted;  f  Set  to  liver  fit  value;  8  Set  to  slowly  per¬ 
fused  tissue  fit  value;  h  Set  to  kidney  fit  value 
#  Fit  value;  *  Value  from  rat  model 

were  euthanized  at  set  time  points  and  plasma,  urine,  liver, 
kidney  and  brain  samples  were  taken. 

Milic  et  al.  (1996)  administered  an  intravenous  dose  of 
16.03  mg  TMB-4/kg  BW  to  male  albino  mice  and  recorded 
the  plasma  time  course.  The  mice  weighed  24-30  g;  a 
midpoint  value  of  27  g  was  used  for  model  simulations. 
This  study  allows  a  limited  evaluation  of  the  rodent 
parameterization  of  the  model  in  which  all  model  inputs 
were  the  same  as  for  the  rat  with  the  exception  of  physi¬ 
ological  parameters.  The  study  also  allows  for  verification 
that  the  model  code  can  simulate  intravenous  dosing. 
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Following  the  rodent  evaluation,  the  model  was 
parameterized  with  two  available  human  data  sets.  Intra¬ 
venous  data  from  Calesnick  and  Christensen  (1967)  were 
useful  in  fitting  the  IM  compartment  PC  and  saturable 
urinary  clearance  parameters  for  humans.  In  this  study,  a 
volunteer  was  administered  24.05  mg  TMB-4/kg  BW 
intravenously.  The  number  of  subjects  was  not  explicitly 
stated,  so  a  single  volunteer  is  assumed  per  time  point,  as 
standard  deviations  are  not  provided  in  the  study.  Male  and 
female  volunteers  ranged  in  age  from  23  to  64  years  old 
with  body  weights  ranging  from  39.5  to  86.2  kg;  a  mid¬ 
point  value  of  62.85  kg  was  used  for  model  simulations. 
The  reported  duration  of  intravenous  infusion  (TInf)  ranged 
from  0.0625  to  0.25  h.  Whole-blood  time  course  samples  at 
various  times  over  an  8-h  period  and  a  24-h  cumulative 
urine  specimen  were  collected. 

A  second  study  was  available  for  use  in  human  param¬ 
eterization.  In  the  Vojvodic  (1970)  study,  male  volunteers 
with  a  mean  body  weight  of  81.89  kg  were  dosed  intra¬ 
muscularly  with  200  mg  of  TMB-4  (approximately 
2.45  mg  TMB-4/kg  BW)  in  2  mL  of  saline.  Whole-blood 
time  course  and  cumulative  urine  amounts  over  4  h  were 
recorded.  This  study  is  useful  in  setting  the  IM  compart¬ 
ment/blood  PC  for  humans  and  fitting  early  time  points  for 
urinary  excretion. 

For  limited  evaluation  of  the  human  parameters,  an 
intramuscular  data  set  from  the  Calesnick  and  Christensen 
(1967)  study  was  utilized.  A  volunteer  was  injected  intra¬ 
muscularly  with  12.02  or  24.05  mg  TMB-4/kg  BW  in 
15  mL  of  sterile  water.  Half  of  the  total  volume  of  the 
injection  (7.5  mL)  was  given  in  each  buttock.  The  length  of 
the  injections  (TIM)  was  set  to  twice  the  value  fit  for  the 
Vojvodic  (1970)  study,  as  it  was  assumed  the  injections 
were  administered  one  after  the  other. 

Monte  Carlo  analyses 

Monte  Carlo  (MC)  analyses  are  used  to  evaluate  variability 
in  model  predictions  and  were  conducted  utilizing  the 
Calesnick  and  Christensen  (1967)  human  intramuscular 
study  parameters.  Physiological  parameters  are  described 
with  a  normal  distribution  while  physicochemical  parame¬ 
ters  (PC,  PA,  VMaxC  and  Km  values)  are  described  with  a 
lognormal  distribution  (Covington  et  al.  2007).  Coefficients 
of  variation  were  set  using  published  values,  except  as  noted 
in  Table  4.  Bounds  were  defined  as  ±  two  standard  devi¬ 
ations  (SD),  except  in  specific  cases.  The  BW  bounds  were 
set  to  the  published  range  for  the  Calesnick  and  Christensen 
(1967)  study  (39.5-86.2  kg)  with  the  analysis  mean  being 
equal  to  the  midpoint  value  of  the  data  range  (62.85  kg). 
Blood  flow  rates  and  tissue  volumes  for  the  slowly  and 
rapidly  perfused  tissues,  as  shown  in  Table  4,  are  not  the 
same  as  those  listed  in  Table  1;  the  mean  for  each  of  these 


values  represents  the  value  calculated  by  the  model  after 
subtracting  the  appropriate  blood  flows  and  tissue  volumes 
from  each  lumped  compartment.  The  bounds  were  then  set 
using  the  post-calculation  number  to  keep  an  MC  analysis 
from  choosing  a  lumped  value  that  is  smaller  than  one  of  the 
component  tissues.  The  metabolic  enzyme  and  urinary 
excretion  affinity  constant  bounds  (Km  and  KmUr, 
respectively)  were  also  set  to  avoid  zero  values,  which 
would  cause  the  Michaelis-Menten  equation  to  divide  by 
zero,  the  saturable  pathways  to  turn  off,  and  metabolism  and 
excretion  to  no  longer  be  concentration  dependent.  The 
model  was  run  for  1,000  iterations  with  all  physiological 
and  physicochemical  parameters  varied  randomly  within 
the  defined  distributions. 

An  MC  analysis  may  select  fractional  values  such  that 
they  no  longer  sum  in  a  physiological  manner  (e.g.,  frac¬ 
tional  blood  flows  summing  to  more  than  1).  Since  mass 
balance  must  be  maintained,  the  randomly  selected  frac¬ 
tional  blood  flows  and  tissue  volumes  are  “rebalanced”  by 
adjusting  each  parameter  in  a  fractional  manner  such  that 
they  sum  to  the  correct  total. 

Sensitivity  analysis 

Sensitivity  analyses  determine  to  which  parameters  the 
model  is  most  responsive.  They  also  help  identify  incon¬ 
sistencies  within  the  model  and  determine  which  critical 
parameters  may  require  further  research  (Kohn  2002; 
Loizou  et  al.  2008).  Sensitivity  analyses  were  conducted  in 
order  to  determine  the  impact  on  selected  endpoints  of 
varying  input  parameters  by  1  %.  These  analyses  were 
conducted  by  increasing  each  parameter  sequentially  by 
1  %  and  calculating  the  resulting  sensitivity  coefficient 
(SC)  for  each  time  point  as  the  increase  in  the  endpoint 
over  the  baseline  endpoint  value  ( i.e .,  with  none  of  the 
input  parameters  varied)  divided  by  the  increase  in 
parameter  value  over  baseline  parameter  value.  Parameters 
were  varied  one  at  a  time  and  were  returned  to  their  ori¬ 
ginal  values  before  the  next  parameter  was  varied.  Dif¬ 
ferences  for  SCs  were  calculated  using  the  central 
difference  method  and  were  normalized  to  both  the 
response  variable  as  well  as  the  parameter  being  varied. 

Sensitivity  analyses  were  run  to  simulate  kinetics  for  3  h 
post-dosing  due  to  the  relatively  short  therapeutic  window 
for  oximes.  Venous  blood  and  brain  concentration  were 
selected  as  the  response  variables  as  sufficient  oxime 
concentrations  in  the  blood  and  target  tissues  are  of  fore¬ 
most  concern  for  achieving  effective  AChE  reactivation. 
Human  study  parameters  from  Calesnick  and  Christensen 
(1967)  and  Vojvodic  (1970)  were  used  to  evaluate  intra¬ 
venous  and  intramuscular  routes,  respectively.  Parameters 
were  considered  to  be  sensitive  if  the  SC  exceeded  0.5,  in 
absolute  value,  meaning  that  the  output  parameter  changed 
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Table  4  Distribution  definitions  for  human  Monte  Carlo  analysis  of  Calesnick  and  Christensen  (1967)  intramuscular  data 

Parameter 

Distribution 

Shape 

Mean 

CV“  (%) 

Bounds 

Lower 

Upper 

Body  weight  (kg) 

BW 

Normal 

62.85b 

30 

39.5° 

86.2° 

Flow  rates  (L/h/kg015) 

Cardiac  output 

QCC 

Normal 

12.89 

9d 

10.57 

15.21 

Ventilation  perfusion  ratio  (QCC/QPC) 

VPR 

Lognormal 

0.962 

15 

0.673 

1.251 

Blood  flow  rates  (fraction  of  cardiac  output) 

Brain 

QBrC 

Normal 

0.114 

20e 

0.068 

0.160 

Diaphragm 

QDC 

Normal 

0.006 

15f 

0.004 

0.008 

Fat 

QFC 

Normal 

0.052 

30 

0.021 

0.083 

Kidneys 

QKC 

Normal 

0.175 

35 

0.053 

0.298 

Liver 

QLC 

Normal 

0.227 

35 

0.068 

0.386 

Rapidly  perfused 

QRC 

Normal 

0.144 

20 

0.086 

0.202 

Slowly  perfused 

QSC 

Normal 

0.282 

15 

0.197 

0.367 

Tissue  volume  (fraction  of  body  weight ) 

Arterial  blood 

VAC 

Normal 

0.0308 

10e 

0.0246 

0.0370 

Venous  blood 

VVC 

Normal 

0.0482 

10e 

0.0386 

0.0578 

Brain 

VBrC 

Normal 

0.020 

10e 

0.016 

0.024 

Diaphragm 

VDC 

Normal 

0.003 

30f 

0.001 

0.005 

Fat 

VFC 

Normal 

0.214 

30 

0.086 

0.342 

Kidneys 

VKC 

Normal 

0.004 

5 

0.0036 

0.0044 

Liver 

VLC 

Normal 

0.026 

5 

0.023 

0.029 

Rapidly  perfused 

VRC 

Normal 

0.013g 

10 

0.010 

0.016 

Slowly  perfused 

VSC 

Normal 

0.627g 

30 

0.251 

0.9° 

Tissue  blood  volume  (fraction  of  tissue  volume) 

Brain 

VBrBC 

Normal 

0.04 

10h 

0.032 

0.048 

Diaphragm 

VDBC 

Normal 

0.01 

30h 

0.004 

0.016 

Fat 

VFBC 

Normal 

0.020 

30h 

0.008 

0.032 

Kidneys 

VKBC 

Normal 

0.36 

5h 

0.32 

0.40 

Liver 

VLBC 

Normal 

0.11 

5h 

0.10 

0.12 

Rapidly  perfused 

VRBC 

Normal 

0.170 

10h 

0.136 

0.204 

Slowly  perfused 

VSBC 

Normal 

0.038 

30h 

0.015 

0.060 

IM  compartment  parameters 

Blood  flow  rate  (L/h/g) 

QIMg 

Normal 

0.0045 

15f 

0.0032 

0.0059 

Tissue  blood  volume  (fraction) 

VIMBC 

Normal 

0.01 

30fh 

0.004 

0.016 

Metabolism  parameters 

Maximum  enzymatic  metabolic  rate  constant  (mg/h/kg0'75) 

VMaxC 

Lognormal 

0.6 

30d 

0.2 

1.0 

Enzyme  affinity  constant  (mg/L) 

Km 

Lognonnal 

0.1 

50d 

0.001° 

0.2 

Urinary  excretion  parameters 

Maximum  saturable  excretion  rate  constant  (mg/h/kg0  75) 

VMaxUrC 

Lognormal 

480.0 

30‘ 

192.0 

768.0 

Excretion  affinity  constant  (mg/L) 

KmUr 

Lognonnal 

900.0 

50‘ 

1.0° 

1800.0 

Tissue/blood  partition  coefficients  (unitless) 

Brain 

PBr 

Lognormal 

0.13 

20 

0.08 

0.18 

Diaphragm 

PD 

Lognormal 

0.94 

20 

0.56 

1.32 

Fat 

PF 

Lognormal 

0.19 

30 

0.08 

0.30 

IM  compartment 

PIM 

Lognonnal 

2.5 

20 

1.5 

3.5 

Kidney 

PK 

Lognormal 

0.94 

20 

0.56 

1.32 

Liver 

PL 

Lognormal 

0.90 

20 

0.5 

1.3 
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Table  4  continued 

Parameter 

Distribution 

Shape 

Mean 

CVa  (%) 

Bounds 

Lower 

Upper 

Rapidly  perfused 

PR 

Lognormal 

0.97 

20 

0.58 

1.36 

Slowly  perfused 

Diffusion  limitation  coefficient  constants  ( IVh/kg °'75) 

PS 

Lognormal 

0.74 

20 

0.44 

1.04 

Brain 

PABrC 

Lognormal 

0.002 

20s 

0.001 

0.003 

Diaphragm 

PADC 

Lognormal 

0.2 

20s 

0.1 

0.3 

Fat 

PAFC 

Lognormal 

0.3 

30s 

0.1 

0.5 

IM  compartment 

PAIMC 

Lognormal 

12.0 

20s 

7.0 

17.0 

Kidney 

PARC 

Lognormal 

1.5 

20s 

0.9 

2.1 

Liver 

PALC 

Lognormal 

0.2 

20s 

0.1 

0.3 

Rapidly  perfused 

PARC 

Lognormal 

1.5 

20s 

0.9 

2.1 

Slowly  perfused 

PASC 

Lognormal 

0.3 

20s 

0.2 

0.4 

a  Coefficient  of  variation  (CV)  values  from  Covington  et  al.  (2007)  unless  otherwise  specified;  b  Calesnick  and  Christensen  (1967)  mid-range 
BW  value;  c  Altered  bounds  (see  text);  d  Clewell  et  al.  (2001);  e  Rapidly  perfused  tissue  CV;  f  Slowly  perfused  tissue  CV;  g  Does  not  represent 
total  flow  or  volume  as  in  Table  1  (see  text);  h  Same  CV  as  for  tissue  volume  compartment;  1  Same  CV  as  for  saturable  metabolism  parameter; 
J  Same  CV  as  for  corresponding  partition  coefficient 


half  as  much  as  the  input  parameter  was  varied.  A 
parameter  was  considered  very  sensitive  if  the  SC  in 
absolute  value  exceeded  1.0. 


Results 

Rat  model  simulations 

The  TMB-4  PBPK  model  was  developed  using  data  from 
Garrigue  et  al.  (1991)  which  were  collected  from  Wistar  rats 
dosed  intramuscularly  with  21 .5  mg  radiolabeled  TMB-4/kg 
BW.  Model  simulations  of  plasma,  tissue  and  cumulative 
urine  time  course  data  are  shown  in  Fig.  3;  overall,  these 
simulations  fit  the  available  data  well  with  minor  exceptions. 
The  model  predicts  a  slightly  earlier  peak  in  the  plasma 
concentration  and  somewhat  over-predicts  the  tail  of  the 
plasma  time  course  (Fig.  3a).  Simulated  brain  and  liver 
concentrations  closely  match  the  available  data  (Fig.  3b,  c) 
while  kidney  concentration  is  over-predicted  at  the  early 
(10  min)  time  point  (Fig.  3b).  Standard  deviations  were  not 
available  for  plasma  or  tissue  data  to  determine  whether  the 
model  predictions  lie  within  the  data  variability.  Predicted 
urinary  elimination  falls  within  one  standard  deviation  of  the 
mean  at  every  time  point,  indicating  reasonable  estimation  of 
urinary  excretion  of  TMB-4  in  rats  (Fig.  3d). 

In  order  to  obtain  these  predictions,  the  model 
required  larger  values  for  kidney/blood,  liver/blood  and 
slowly  perfused  tissue  compartment/blood  PCs  (100.0, 
12.0  and  10.0,  respectively)  versus  the  smaller  predicted 
values  (0.84,  0.69  and  0.71,  respectively).  Conversely, 


the  model  required  a  smaller  brain/blood  PC  (0.13) 
compared  to  the  predicted  value  (0.91).  The  length  of  the 
injection  (TIM),  the  PA  values  and  the  tissue/blood  PC 
for  the  IM  compartment  were  also  fit  for  the  Garrigue 
et  al.  (1991)  data. 

Evaluation  of  model  and  rodent  parameters 

Since  only  two  sets  of  rodent  time  course  kinetic  data  for 
TMB-4  were  found  in  the  literature,  the  remaining  data  set 
(Milic  et  al.  1996)  was  used  for  limited  evaluation  of  the 
model.  Milic  et  al.  (1996)  published  plasma  time  course 
concentrations  for  mice  dosed  intravenously  with  16.03  mg 
TMB-4/kg  BW.  The  model  is  able  to  predict  these  plasma 
time  course  data  (Fig.  4)  utilizing  mouse  physiological 
parameters  (Table  2)  and  rat  physicochemical  parameters 
(Table  3).  It  was  necessary  to  fit  the  length  of  intravenous 
infusion  (TInf),  as  this  is  not  stated  within  Milic  et  al. 
(1996).  The  injection  duration  was  set  to  0.01  h.  Over¬ 
prediction  of  the  first  time  points  occurs  because  the  model 
assumes  instant  mixing  of  venous  blood.  The  model  under¬ 
predicts  mouse  plasma  TMB-4  levels  in  the  tail  of  the  time 
course.  Overall,  the  evaluation  of  the  model  indicates  a 
successful  prediction  utilizing  the  rat  parameters  for  a 
similar  species  (mouse). 

Human  model  simulations 

The  model  demonstrates  its  ability  to  be  used  for  interspe¬ 
cies  extrapolation  by  predicting  whole-blood  concentrations 
and  cumulative  urinary  elimination  data  for  humans  with 
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Fig.  3  Model  simulations 
(lines)  for  a  plasma,  b  kidney 
(dashed  line  and  open 
triangles ),  liver  (solid  line  and 
closed  triangle)  and  c  brain 
concentrations,  and 
d  cumulative  urinary  excretion 
compared  with  measured  rat 
data  (symbols,  mean  ±  1  SD) 
from  Garrigue  et  al.  (1991). 
Rats  were  administered  a  single 
intramuscular  dose  of  Re¬ 
labeled  TMB-4  equivalent  to 
21.5  mg/kg  BW 


Fig.  4  Model  simulation  (line)  of  plasma  concentration  compared 
with  measured  mouse  data  (symbols,  mean  ±  1  SD)  from  Milic  et  al. 
(1996).  Mice  received  an  intravenous  dose  of  16.03  mg  TMB-4/kg 
BW 

relatively  few  parameter  adjustments.  Intravenous  data 
from  Calesnick  and  Christensen  (1967)  and  intramuscular 
data  from  Vojvodic  (1970)  were  used  to  parameterize  the 
model  for  human  kinetics.  The  change  in  species  necessi¬ 
tated  alterations  to  the  maximum  rate  constants  for  metab¬ 
olism  and  excretion  (VMaxC,  VMaxUrC).  In  addition,  dose 
route  parameters,  including  the  IM  compartment/blood  PC 
(PIM)  and  the  lengths  of  intramuscular  and  intravenous 


infusions  (TIM,  TInf),  are  altered  to  fit  not  only  the  changed 
physiology,  but  also  study  protocol  differences. 

Predicted  PC  values  were  largely  utilized  for  the  human 
parameterization  of  the  TMB-4  model.  Tissue/blood  PC 
values  for  the  liver,  kidney  and  slowly  perfused  tissues, 
previously  fit  to  the  rat  data,  were  returned  to  QSPR-pre- 
dicted  levels  for  the  human  data  sets.  The  brain/blood  PC 
remained  the  same  as  the  fit  value  (0.13)  from  the  Garrigue 
et  al.  (1991)  rat  data  due  to  the  lack  of  human  brain  data. 
The  IM  compartment/blood  PC  was  fit  using  the  Vojvodic 
(1970)  data;  visual  fitting  to  these  intramuscular  data 
results  in  a  tissue/blood  PC  of  2.5  for  the  IM  compartment. 
The  human  parameterization  used  rat  PA  values. 

The  model  is  able  to  simulate  the  Calesnick  and  Chris¬ 
tensen  (1967)  intravenous  whole-blood  data  well  (Fig.  5a). 
The  volunteer  was  administered  24.05  mg  TMB-4/kg  BW 
intravenously.  The  length  of  injection  (TInf)  is  set  to  the 
minimum  study  value  (0.0625  h)  as  it  best  simulated  the  data 
curve.  The  model  under-predicts  the  mean  24-h  urinary  output 
by  12  %  (Fig.  5b).  Rather  than  fitting  excretion  parameters  to 
this  single  data  point,  the  data  from  Vojvodic  (1970)  were 
utilized  to  set  the  parameters  for  this  endpoint. 

The  Vojvodic  (1970)  data  set  was  utilized  for  fitting 
human  excretion  rates  and  intramuscular  administration 
parameter  values.  The  model  fits  the  whole-blood  time 
course  and  the  cumulative  urinary  elimination  data  well 
(Fig.  6).  Standard  deviation  values  were  only  available  for 
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the  blood  time  course  and  illustrate  the  variability  of  TMB-4 
blood  concentrations  between  individuals.  The  model  is 
capable  of  simulating  the  urinary  excretion  curve  at  early 
time  points  which  cannot  be  determined  for  rats  due  to  the 
lack  of  early  data. 

Evaluation  of  human  parameters  and  Monte  Carlo 
analyses 

To  evaluate  the  human  TMB-4  model  parameterization,  an 
intramuscular  data  set  of  TMB-4  blood  levels  from  the 
Calesnick  and  Christensen  (1967)  study  was  reserved.  In  the 
intramuscular  portion  of  this  study,  a  volunteer  was 
administered  12.02  or  24.05  mg  TMB-4/kg  BW.  The  model 
reasonably  predicts  the  low  dose  whole-blood  data  (Fig.  7a) 
with  the  exception  of  the  curious  flatness  of  the  data  at  early 
time  points  not  seen  previously  in  the  rat  (Fig.  3)  or  the 
human  (Fig.  6)  intramuscular  data  of  Vojvodic  (1970).  The 
model  uniformly  under-predicts  the  higher  intramuscular 
dose  by  approximately  20  %  (Fig.  7b). 

To  evaluate  the  effect  of  normal  human  variability  on 
model  predictions,  Monte  Carlo  analyses  were  performed. 
Calesnick  and  Christensen  (1967)  accepted  1 1  male  and  29 
female  volunteers  ranging  greatly  in  age  (23-64  years)  and 
body  weight  (39.5-86.2  kg).  Figure  7  shows  the  results  of 
the  MC  analyses  of  the  Calesnick  and  Christensen  (1967) 
human  intramuscular  study  using  the  study  BW  range 
(39.5-86.2  kg)  as  bounds  for  that  parameter  distribution 
and  the  midpoint  of  that  range  (62.85  kg)  as  the  mean  for 
the  distribution.  The  MC  estimates  show  that  the  observed 
data  from  one  or  more  volunteers  per  dose  fell  within  2  SD 
of  model  predicted  values,  reflecting  a  fit  within  95  %  of 
normal  human  population  variability  given  the  body  weight 
range  specified. 

Calesnick  and  Christensen  (1967)  report  a  24-h  elimi¬ 
nation  of  just  34  %  for  both  the  high  and  low  doses.  This 
low  level  of  elimination  is  not  seen  in  the  rat  data  above, 
nor  in  the  human  intravenous  and  other  human 


intramuscular  data  sets  previously  presented;  therefore, 
predictions  of  these  apparently  anomalous  data  are  not 
attempted. 

Sensitivity  analysis 

The  model  is  sensitive  to  several  parameters  over  a  3-h 
simulation.  The  endpoints  of  concern  for  the  sensitivity 
analysis  are  TMB-4  blood  and  brain  concentrations.  Both 
intravenous  and  intramuscular  dosing  routes  were  ana¬ 
lyzed;  all  physiological  and  physicochemical  inputs  into 
the  model  were  varied.  SC  ranges  are  listed  for  each  sen¬ 
sitive  and  very  sensitive  parameter  by  route  and  endpoint  in 
Table  5. 

The  blood  concentration  of  TMB-4  administered  by  the 
intravenous  route  is  sensitive  to  six  parameters.  Figure  8a 
shows  two  of  the  physiological  parameters  to  which  the 
blood  concentration  is  sensitive:  the  volume  of  the  slowly 
perfused  tissue  compartment  (VSC)  and  the  volume  of 
venous  blood  (VVC).  The  model  is  sensitive  to  VVC  only 
during  the  initial  distribution  phase;  the  model  is  sensitive 
to  VSC  initially  and  again  during  the  end  of  the  critical  3-h 
time  period.  Figure  8b  shows  the  SC  time  courses  of  an 
administration  parameter,  length  of  intravenous  injection 
(TInf),  to  which  the  model  is  sensitive  only  during  the 
administration  of  TMB-4.  Three  physicochemical  param¬ 
eter  SC  time  courses  (the  kidney/blood  PC  (PK),  the 
maximum  saturable  excretion  rate  constant  (VMaxUrC) 
and  the  excretion  affinity  constant  (KmUr))  are  also  shown 
in  Fig.  8b.  These  parameters  affect  kidney  uptake  and 
excretion,  and  their  variance  results  in  model  sensitivity  in 
the  latter  half  of  the  3-h  time  course. 

The  brain  concentration  of  TMB-4  by  the  intravenous 
route  is  sensitive  to  the  same  physiological  and  adminis¬ 
tration  parameters  as  the  blood  concentration.  Again,  this 
endpoint  is  sensitive  to  VVC  (Fig.  9a)  and  TInf  (Fig.  9b) 
during  the  initial  distribution  of  TMB-4;  however,  brain 
concentration  is  sensitive  to  VSC  for  almost  the  entire  3-h 


Fig.  5  Model  simulations 
(lines)  compared  with  measured 
human  intravenous  data 
(symbols)  for  a  whole-blood 
time  course  and  b  cumulative 
urinary  excretion  from 
Calesnick  and  Christensen 
(1967).  Volunteers  were 
administered  24.05  mg  TMB-4/ 
kg  BW 
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Fig.  6  Model  simulations  (lines)  compared  with  measured  human 
intramuscular  data  (symbols,  mean  ±  1  SD,  when  available)  for 
a  whole-blood  time  course  and  b  cumulative  urinary  excretion  from 


A 


Fig.  7  Monte  Carlo  estimates  of  variability  of  blood  concentration 
(lines)  compared  with  measured  human  intramuscular  data  (symbols) 
from  Calesnick  and  Christensen  (1967).  Average  predictions  from 
1000  simulations  (solid  line)  are  shown  as  well  as  ±  two  SD 

time  period  (Fig.  9a).  Brain  concentration  was  also  sensi¬ 
tive  to  five  additional  parameters.  Sensitivity  is  shown  to 
three  more  physiological  parameters  (cardiac  output 
(QCC),  volume  of  arterial  blood  (VAC)  and  volume  of  the 
brain  (VBrC))  (Fig.  9a).  Sensitivity  to  QCC  and  VAC 
occur  only  initially,  while  brain  concentration  is  sensitive 
to  VBrC  throughout  most  of  the  3  h.  Variance  of  two 
physicochemical  parameters  (the  brain/blood  PC  (PBr)  and 
brain  diffusion  limitation  coefficient  constant  (PABrC)) 
result  in  sensitivity  of  this  endpoint  as  well  (Fig.  9b).  The 
model  is  sensitive  to  PABrC  throughout  most  of  the 
3  hours,  while  it  is  sensitive  to  PBr  only  at  the  very  end  of 
the  evaluated  time  period.  Brain  concentration  output  was 
found  to  be  very  sensitive  to  only  two  of  these  eight 
parameters,  VSC  and  VBrC,  by  the  criteria  of  SC  values 
greater  in  absolute  value  than  1.0. 


B 


Vojvodic  (1970).  Male  volunteers  were  dosed  intramuscularly  with 
2.45  mg  TMB-4/kg  BW 


B 


predictions  (dotted  line)  which  encompass  95  %  of  human  variability. 
Volunteers  were  administered  intramuscular  doses  of  a  12.02  or 
b  24.05  mg  TMB-4/kg  BW 

The  blood  concentration  endpoint  for  the  intramuscular 
route  was  found  to  be  sensitive  to  only  four  parameters. 
Again,  VVC  and  VSC  were  two  of  the  physiological 
parameters;  this  endpoint  was  sensitive  to  VVC  at  the  very 
beginning  of  distribution  of  the  drug  and  was  sensitive  to 
VSC  during  the  first  and  last  portions  of  the  time  course 
(Fig.  10).  Blood  concentration  was  sensitive  to  another 
physiological  parameter  (blood  flow  to  the  IM  compart¬ 
ment  (QIMg))  and  a  physicochemical  parameter  (the  IM 
compartment/blood  PC);  variance  of  these  two  route-spe¬ 
cific  parameters  resulted  in  sensitivity  only  during  the  first 
hour  of  the  time  course.  The  model  was  very  sensitive  to 
variance  of  only  one  parameter,  VSC. 

Finally,  the  brain  concentration  of  TMB-4  following 
intramuscular  administration  was  sensitive  to  12  parame¬ 
ters.  As  with  all  of  the  other  analyses,  brain  concentration 


DISTRIBUTION  A.  Approved  for  public 


<£)  Springer 

release;  distribution  unlimited. 


672 


Arch  Toxicol  (2013)  87:661-680 


Table  5  Minimum  and  maximum  sensitivity  coefficient  (SC)  values  for  human  intravenous  and  intramuscular  sensitivity  analyses 
Parameter  Endpoint 


Intravenous  concentration  Intramuscular  concentration 


Blood 

Brain 

Blood 

Brain 

Body  weight 

BW 

NS 

NS 

NS 

-0.624,  0.060 

Flow  rates 

Cardiac  output 

QCC 

NS 

-0.085,  0.889 

NS 

-0.082,  1.52 

Blood  flow  rates 

Blood  flow  to  brain 

QBrC 

NS 

NS 

NS 

0,  0.671 

Blood  flow  to  IM  compartment 

QIMg 

NA 

NA 

-0.413,  0.999 

0.255,  0.999 

Blood  flow  to  rapidly  perfused  tissues 

QRC 

NS 

NS 

NS 

-0.044,  0.554 

Tissue  volume 

Volume  of  arterial  blood 

VAC 

NS 

-0.633,  0.008 

NS 

-0.756,  0.009 

Volume  of  brain 

VBrC 

NS 

-1.32,  -0.435 

NS 

-1.65,  -0.585 

Volume  of  slowly  perfused  tissues 

VSC 

-0.579,  0.620 

0.411,  2.49 

-0.590,  1.62 

0.821,  3.92 

Tissue  blood  volume 

Volume  of  venous  blood 

VVC 

-0.620,  0.122 

-0.668,  0.012 

-0.786,  0.063 

-0.710,  0.013 

Volume  of  brain  blood 

VBrBC 

NS 

NS 

NS 

-0.625,  0.032 

Urinary  excretion  parameters 

Maximum  saturable  excretion  rate  constant 

VMaxUrC 

-0.771,  0 

NS 

NS 

NS 

Excretion  affinity  constant 

KmUr 

0,  0.752 

NS 

NS 

NS 

Tissue/blood  partition  coefficients 

Brain/blood  PC 

PBr 

NS 

0,  0.551 

NS 

NS 

IM  compartment/blood  PC 

PIM 

NA 

NA 

-0.987,  0.414 

-0.987,  -0.254 

Kidney/blood  PC 

PK 

-0.744,  -0.001 

NS 

NS 

NS 

Diffusion  limitation  coefficient  constants 

Brain  diffusion  limitation 

Coefficient  constant 

PABrC 

NS 

0.444,  0.998 

NS 

0.604,  0.999 

Administration  route  parameter 

Duration  of  injection  (intravenous  or  intramuscular) 

TInf  or  TIM 

-0.991,  0.505 

-0.990,  0.983 

NS 

NS 

Bold  SC  values  are  >  10.51;  Underlined  SC  values  are  >  11.01 
NA  not  applicable  to  the  route,  NS  not  sensitive 
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Fig.  8  Time  courses  of  sensitivity  coefficients  for  six  model  parameters  to  which  model  predictions  of  TMB-4  blood  concentration  are  sensitive 
following  intravenous  administration 
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Fig.  9  Time  courses  of  sensitivity  coefficients  for  eight  model  parameters  to  which  model  predictions  of  TMB-4  brain  concentration  are 
sensitive  following  intravenous  administration 


Time  (hours) 

Fig.  10  Time  courses  of  sensitivity  coefficients  for  four  model 
parameters  to  which  model  predictions  of  TMB-4  blood  concentration 
are  sensitive  following  intramuscular  administration 

was  sensitive  to  changes  in  VVC  and  VSC  (Fig.  11a).  This 
endpoint  was  sensitive  to  VVC  early  in  the  time  course  and 
was  sensitive  to  VSC  throughout  the  3-h  period.  The 
endpoint  was  sensitive  to  QCC  and  VAC  (Fig.  11a)  in  the 
manner  previously  seen  in  Fig.  9a  for  the  intravenous 
analysis.  Sensitivity  was  also  shown  to  QIMg  and  PIM,  in 
an  equal  and  opposite  manner  (Fig.  1  lb),  as  shown  for  the 
blood  endpoint  in  Fig.  10.  The  intramuscular  endpoint 
analysis  revealed  sensitivity  to  two  additional  physiologi¬ 
cal  parameters:  the  blood  flow  to  the  rapidly  perfused  tissue 
compartment  (QRC)  and  BW  (Fig.  lib).  Both  were  only 
significant  in  the  early  distribution  stage  of  the  3-h  analy¬ 
sis.  Finally,  the  brain  concentration  of  TMB-4  was  sensi¬ 
tive  to  four  brain  parameters  (Fig.  11c),  of  which  only 
PABrC  was  seen  in  the  intravenous  analysis  (Fig.  9b).  The 
intramuscular  analysis  showed  sensitivity  to  blood  flow  to 
the  brain  (QBrC)  and  volume  of  blood  in  brain  tissue 
(VBrBC),  both  early  in  the  time  course  (Fig.  lie).  The 


volume  of  the  brain  itself  (VBrC)  was  a  sensitive  parameter 
throughout  the  time  course  in  this  analysis.  Variance  of 
three  parameters,  QCC,  VBrC  and  VSC,  resulted  in  very 
sensitive  SC  values  by  the  chosen  criterion. 

Discussion 

The  model  presented  here  is  the  first  published  PBPK 
model  of  an  oxime.  Previous  kinetic  models  of  oximes 
utilized  one-  and  two-compartment  models  to  simulate 
blood  concentrations  for  individual  studies  (Stemler  et  al. 
1991;  Worek  et  al.  2005).  A  PBPK  model  has  the  advan¬ 
tage  of  interspecies  and  cross-route  extrapolation.  This 
PBPK  model  was  initially  developed  to  predict  plasma  and 
tissue  time  course  and  urinary  elimination  data  following 
intramuscular  administration  of  the  oxime  TMB-4  in  rats. 
The  model  was  then  evaluated  with  mouse  intravenous 
data.  Following  initial  development,  the  model  was 
parameterized  for  humans  utilizing  intravenous  and  intra¬ 
muscular  data  and  evaluated  using  additional  human 
intramuscular  data.  The  PBPK  model  structure  incorpo¬ 
rates  multiple  tissues  important  to  the  administration,  dis¬ 
tribution,  metabolism  and  elimination  of  oximes.  The 
structure  also  includes  target  tissues  for  OP  action  and 
oxime  reactivation  of  AChE  activity  (brain,  diaphragm), 
which  will  be  useful  in  tying  the  TMB-4  model  to  an  OP 
model  for  future  PBPK  and  pharmacodynamic  model 
development. 

Also  incorporated  in  the  model  structure  is  an  intra¬ 
muscular  dosing  compartment.  The  IM  compartment  is 
different  from  the  other  model  compartments  in  two  ways: 
(1)  it  is  created  within  the  model  code  only  with  an  intra¬ 
muscular  dose  and  (2)  the  model  calculates  the  size  of  the 
compartment  based  on  the  volume  of  the  dose.  The  design 
of  this  compartment  allows  the  model  to  more  correctly 
simulate  different  dose  volumes,  which  ranged  in  humans 
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Fig.  11  Time  courses  of  sensitivity  coefficients  for  12  model  parameters  to  which  model  predictions  of  TMB-4  brain  concentration  are  sensitive 
following  intramuscular  administration 


from  2  mL  (Vojvodic  1970)  to  15  mL  (Calesnick  and 
Christensen  1967).  Dosing  and  absorption  via  the  IM 
compartment  are  dependent  upon  the  length  of  time  for  the 
injection  (TIM),  the  volume  of  the  injection  (VolIM),  the 
blood  flow  to  the  injection  site  (QIMg)  and  the  partition 
coefficient  of  the  IM  compartment  (PIM).  Separate  values 
for  PIM  were  fit  for  rodents  and  humans.  A  single  human 
value  was  used  to  simulate  the  model  development  data 
(Vojvodic  1970)  and  predict  the  limited  evaluation  data 
(Calesnick  and  Christensen  1967  data)  with  reasonable 
results. 

The  IM  compartment  designed  for  this  model  is  unlike 
those  in  previously  published  models.  The  DFP  model  on 
which  this  project  was  based  (Gearhart  et  al.  1994)  added 
the  dose  directly  into  the  slowly  perfused  tissue  compart¬ 
ment,  using  an  artificially  slow  injection  rate  to  regulate  the 
speed  of  chemical  absorption.  Another  model  utilized  a 
sub-compartment  of  the  slowly  perfused  tissues  for  the 
purpose  of  intramuscular  dosing  (Chen  and  Seng  2012); 
however,  it  modeled  absorption  from  this  sub-compartment 
as  going  to  the  larger  slowly  perfused  tissues  compartment 
rather  than  directly  into  the  blood,  which  would  still 
require  the  drug  concentration  to  come  to  equilibrium  with 
the  entire  slowly  perfused  tissues  compartment  before 
partitioning  into  the  blood  stream.  Both  of  these  methods 


are  less  physiologically  accurate  as  this  compartment  is  a 
lumped  entity  comprised  of  various  tissues  throughout  the 
body  rather  than  in  a  single  location  ( e.g .,  liver);  therefore, 
in  the  model  described  herein,  the  compound  enters  the 
bloodstream  prior  to  reaching  equilibrium  in  this  diffuse 
tissue. 

Craigmill  (2003)  and  Yang  et  al.  (2012)  utilized  similar 
stand-alone  injection  site  compartments  for  intramuscular 
administration  models  in  sheep  and  pigs,  respectively.  The 
sizes  of  these  compartments  were  fixed  at  a  set  volume  (0.5 
or  1  kg,  respectively)  over  which  the  volume  of  an  injec¬ 
tion  is  equilibrated  prior  to  diffusion  into  the  blood  com¬ 
partment.  With  the  varied  dose  sizes  used  in  the  TMB-4 
studies  simulated  here,  utilizing  a  fixed  volume  method 
still  has  the  dose  being  injected  into  a  relatively  large 
compartment,  and  results  in  a  diluted  dose  absorbing  into 
the  blood  stream.  This  could  require  re-adjusting  parame¬ 
ters  governing  uptake  into  the  blood  whenever  the  size  of 
the  dose  changes  significantly.  By  using  a  volume  for  the 
IM  compartment  that  is  determined  by  the  size  of  the  dose, 
the  concentration  being  absorbed  into  the  blood  stream  will 
be  more  consistent  for  varying  doses  and  should  require 
fewer  adjustments  to  the  uptake  parameters.  Thus,  the  IM 
compartment  used  in  this  PBPK  model  has  the  advantage 
of  adapting  to  real-world  administration  routes  (the  IM 
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compartment  is  created  only  when  the  intramuscular  route 
is  utilized),  injection  rates,  and  dosing  volumes,  potentially 
without  having  to  adjust  the  fit  PIM  value  between  same- 
species  studies. 

Metabolism  is  included  in  model,  contrary  to  the  conclu¬ 
sions  reached  by  Garrigue  et  al.  (1991),  the  study  on  which  this 
model  was  built.  Garrigue  and  colleagues  utilized  chroma¬ 
tography  to  determine  that  no  metabolites  of  TMB-4  were 
present  in  rat  urine;  however,  their  data  indicated  incomplete 
excretion  of  TMB-4  at  72  h,  long  after  the  oxime  had  left  the 
blood.  Metabolism  of  TMB-4  had  been  documented  in  rats, 
with  de  Miranda  et  al.  (1972)  first  identifying  the  minor 
metabolite,  trimethylene- 1  -(4-aldoximinopyridinium)- 1  ’  -(4- 
cyanopyridinium)  ion  (TACN)  and  Morgan  et  al.  (1982) 
identifying  the  metabolite,  trimethylene-  l-(4-aldoximinopy- 
ridinium)-r-(4-carbox-amidopyridinium)  ion  (TACARB)  in 
rat  urine.  For  the  Garrigue  et  al.  ( 1 99 1 )  data,  the  TMB-4  PBPK 
model  could  not  replicate  the  blood,  tissue  and  urine  data 
without  metabolism  to  eliminate  some  of  the  TMB-4.  The 
model  simulations  were  visually  fit  to  the  Garrigue  data  set 
using  an  estimated  15  %  metabolism  over  72  h. 

The  model  structure  also  incorporates  saturable  urinary 
clearance  using  the  Michaelis-Menten  equation  to  describe 
TMB-4  elimination.  First-order  urinary  clearance  does  not 
adequately  describe  the  rapid  excretion  of  TMB-4,  sug¬ 
gesting  a  potentially  different  mechanism  of  uptake  in 
addition  to  absorption  in  the  renal  tubules  (Sidell  and  Groff 
1971;  Swartz  and  Sidell  1974).  Kayouka  et  al.  (2011)  have 
found  that  decreasing  organic  cation  transport  in  rats, 
through  competitive  pretreatment  with  the  organic  cation 
tetraethylammonium  or  in  knockout  mice  deficient  in 
OCTs  1  and  2,  resulted  in  significantly  slower  excretion  of 
the  oxime  2-PAM.  Located  primarily  in  proximal  tubule 
cells,  OCT  transport  from  the  blood  across  the  basolateral 
membrane  is  dependent  on  a  negative  intracellular  elec¬ 
trical  gradient  (Shu  2011),  thereby  rendering  this  transport 
saturable  and  subject  to  competitive  interaction  from  other 
cationic  drugs  (Choi  and  Song  2008;  Zhang  et  al.  1998). 

Diffusion-limited  uptake  is  utilized  in  all  tissues  and 
requires  a  tissue/blood  PC  and  diffusion  limitation  coeffi¬ 
cient  constants  (PA  values)  for  each  compartment.  Fitting 
the  data  from  Garrigue  et  al.  (1991)  requires  larger  tissue/ 
blood  PCs  for  kidney,  liver  and  the  slowly  perfused  tissue 
compartment  (100.0,  12.0  and  10.0,  respectively)  than 
those  predicted  by  QSPR,  which  predicted  uptake  to  be 
similar  across  these  tissues  (Table  1),  based  on  the  organ 
lipid  and  water  contents.  Preferential  uptake  of  oximes  has 
been  seen  in  rat  liver  and  kidney  by  previous  investigators 
(Lundy  et  al.  1990).  Differences  between  predicted  and  fit 
PC  values  are  thought  to  be  the  result  of  PC  predictions  not 
accounting  for  uptake  limitation  by  active  transport  efflux 
mechanisms  (Schmitt  2008),  in  addition  to  possible  effects 
of  tissue  and/or  blood  binding  (see  Species  Differences 
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below).  The  model  is  also  unable  to  simulate  the  data  with 
the  predicted  brain/blood  PC  of  0.91  in  rats.  Decreasing  the 
PC  to  0.13  is  not  unreasonable  due  to  transporter-limited 
diffusion  across  the  blood-brain  barrier  (BBB).  Oximes  are 
substrates  for  the  active  P-glycoprotein  efflux  transporter 
located  in  the  BBB  and  do  not  easily  enter  the  brain 
(Joosen  et  al.  2011).  The  model  simulation  of  the  Garrigue 
et  al.  (1991)  study  indicates  that  the  maximum  concentra¬ 
tion  in  the  plasma  is  nearly  18  times  that  in  the  brain 
(Fig.  3a,  c);  therefore,  the  necessity  of  the  lower  blood/ 
brain  PC  was  evident  during  model  parameterization  in 
order  to  simply  simulate  the  efflux  of  the  oxime  through  the 
BBB. 

The  model  also  requires  small  PA  values,  indicative  of 
slow  entry  of  TMB-4  into  tissues.  Although  very  water 
soluble,  the  charge  and  large  size  of  TMB-4  molecules  could 
slow  their  movement  into  tissues.  The  brain  compartment 
requires  the  smallest  PA  value,  necessary  to  simulate  the 
influence  of  the  blood-brain  barrier  on  TMB-4  kinetics. 

Species  differences 

The  TMB-4  PBPK  model  proves  useful  in  highlighting 
differences  in  species’  kinetics  for  this  and  potentially 
other  oximes.  The  rat  parameterization  of  this  model 
requires  kidney/blood,  liver/blood  and  slowly  perfused 
tissue/blood  PCs  much  higher  than  those  predicted  by 
QSPR,  while  the  brain/blood  PC  needs  to  be  much  lower. 
Without  the  fit  PCs,  the  model  cannot  attain  the  relatively 
high  concentrations  of  TMB-4  in  the  liver  and  kidney  nor 
the  very  low  concentrations  of  TMB-4  in  the  brain  seen  in 
data  from  Garrigue  et  al.  (1991).  The  ability  of  the  model 
to  reasonably  predict  mouse  data  from  Milic  et  al.  (1996) 
while  only  fitting  the  length  of  intravenous  injection  (TInf) 
indicates  that  the  kinetics  in  rodents  are  potentially  similar. 
Human  data,  however,  cannot  be  simulated  utilizing  these 
adjusted  PCs.  For  the  human  parameterization,  kidney/ 
blood,  liver/blood  and  slowly  perfused  tissue/blood  PCs  are 
adjusted  back  down  to  the  QSPR-predicted  values.  This 
change  accomplished  satisfactory  simulations  and  predic¬ 
tions  to  human  blood  and  urine  data,  as  human  tissue  data 
are  not  attainable.  Potential  sources  of  what  may  be  sig¬ 
nificant  kinetic  differences  between  rodents  and  humans 
are  explored  below. 

Garrigue  et  al.  (1991)  remarked  on  the  high  affinity  of 
the  oximes  they  tested,  including  TMB-4,  for  kidney,  liver, 
cartilage  and  intervertebral  discs  in  the  rat.  A  year  later, 
Garrigue  and  other  colleagues  (Maurizis  et  al.  1992)  pub¬ 
lished  an  in  vitro  binding  study  of  TMB-4  with  rabbit 
cartilaginous  tissue  cultures.  This  study  resulted  in  TMB-4 
ionic  binding  not  with  the  chondrocytes  and  fibroblasts 
themselves,  but  with  the  proteoglycans  produced  by  these 
structures.  Thus,  the  Maurizis  et  al.  (1992)  study  may 
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explain,  in  part,  the  necessity  of  assuming  an  increased 
TMB-4  affinity  for  slowly  perfused  tissues,  which  include 
cartilage  and  similar  structures,  expressed  in  the  model  as  a 
high  tissue/blood  PC.  No  data  were  available  for  deter¬ 
mining  whether  human  cartilage  exhibits  these  same 
kinetics. 

The  affinity  of  TMB-4  for  rat  kidney  and  liver  seen  in 
the  Garrigue  et  al.  (1991)  study  may  be  due,  at  least  in  part, 
to  the  affinity  of  oximes  for  OCTs,  which  are  expressed  in 
both  liver  and  kidney  (Dresser  et  al.  2000)  with  kinetics 
and  selectivity  for  particular  organic  cations  differing 
between  species.  Dresser  et  al.  (2000)  especially  noted 
differences  in  organic  cation  uptake  affinity  and  kinetics 
between  Xenopus  laevis  oocytes  expressing  human  OCT  1 
versus  oocytes  expressing  rat  or  mouse  OCT  1.  Species- 
specific  affinity  for  this  and  other  OCTs  could  contribute  to 
the  need  for  significantly  different  liver  and  kidney  tissue/ 
blood  PCs  when  switching  from  rodents  to  humans. 

A  literature  search  for  blood-binding  information  with 
TMB-4  was  unproductive.  A  search  for  blood  binding  with 
MMB-4  also  failed  to  net  articles.  Binding  research  results 
for  other  oximes  indicate  potential  species  differences. 
2-PAM  has  been  shown,  in  vitro,  to  be  bound  49  to  67  % 
with  sheep  plasma  proteins  (Srivastava  and  Malik  1988) 
and  64.7  %  with  goat  plasma  proteins  (Rahal  and  Malik 
2011);  however,  2-PAM  binding  was  not  observed  in  vitro 
in  human  plasma  or  rat  brain  homogenate  (Firemark  et  al. 
1964).  Ionic  binding  with  a  negatively  charged  structure  is 
theoretically  possible  for  any  oxime,  leaving  potential 
binding  as  a  research  gap  for  individual  oximes. 

No  attempt  was  made  to  alter  the  brain/blood  PC  in 
humans  from  the  small  value  which  had  been  fit  to  avail¬ 
able  rat  brain  data,  although  the  other  fit  PCs  from  the  rat 
model  were  not  utilized  in  the  human  parameterization 
(i.e.,  liver,  kidney  and  slowly  perfused  tissues/blood  PCs) 
due  to  assumed  species  differences.  The  homology  of  the 
rat  and  human  blood-brain  barrier  to  TMB-4  is  not  spe¬ 
cifically  known;  however,  drug  efficacy  limitations  due  to 
P-glycoprotein  and  other  efflux  transporters  in  the  BBB  are 
well  studied  in  both  humans  and  rats.  Blocking  these 
membrane  efflux  transporters  has  been  shown  to  improve 
brain  drug  concentrations  in  studies  utilizing  rats  as  the 
animal  model  (Joosen  et  al.  2011;  Miller  et  al.  2008).  In  the 
absence  of  human  brain  data,  the  assumption  was  made 
that  the  fit  value  for  the  rat  blood/brain  PC  is  similar 
enough  to  the  human  blood/brain  PC.  Overall,  model  pre¬ 
dictions  for  three  species  indicate  a  structure  capable  of 
interspecies  extrapolation  of  TMB-4  kinetics. 

Monte  Carlo  analyses 

Following  initial  verification  of  the  human  model  param¬ 
eterization,  MC  analyses  were  conducted  due  to  the  poor 
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prediction  for  one  of  the  two  human  evaluation  data  sets. 
The  model  is  successful  in  predicting  the  lower  dose  of  the 
intramuscular  data  from  Calesnick  and  Christensen  (1967), 
but  significantly  under-predicts  the  higher  dose.  It  is 
unclear  from  the  methods  published  in  Calesnick  and 
Christensen  (1967)  whether  more  than  one  volunteer  was 
used  at  each  dose  level.  It  is  entirely  possible  that  the 
volunteer  used  for  the  upper  dose  was  physiologically 
different  than  the  lower-dose  individual  or  the  mid-range 
bodyweight  and  physiological  parameters  utilized  for  the 
model  predictions.  MC  analyses  were  used  to  look  at  var¬ 
iability  of  all  physiological  parameters,  since  the  use  of 
male  and  female  volunteers  with  the  range  in  age  and 
bodyweight  described  in  the  study  indicate  a  physiologi¬ 
cally  diverse  group.  Following  the  MC  analyses,  the  high- 
dose  data  are  shown  to  fall  within  two  SD  of  the  human 
variability  likely  seen  within  this  study.  Use  of  MC  anal¬ 
yses  proved  insightful,  demonstrating  that  variation  in 
physiological  and  physiochemical  parameters  between 
individuals  impacts  the  model  predictions  and  is  found  to 
be  consistent  with  the  TMB-4  blood  clearances  seen  in  the 
Calesnick  and  Christensen  (1967)  high-dose  data. 

Sensitivity  analyses 

As  a  final  step  in  building  this  model,  sensitivity  analyses 
were  run.  Sensitivity  analyses  help  characterize  how  well  the 
model  reflects  the  system  in  question;  the  parameters  to 
which  the  model  is  sensitive  should  make  sense  physiologi¬ 
cally  and  kinetically  (Loizou  et  al.  2008).  The  analyses  begin 
with  the  choice  of  critical  endpoints,  in  this  case  TMB-4 
concentrations  in  blood  and  brain.  TMB-4  is  eliminated  very 
quickly,  which  reduces  its  efficacy  in  reactivating 
OP-inhibited  AChE  (Kayouka  et  al.  201 1).  A  3-h  time  course 
was  chosen  to  capture  the  portion  of  the  concentration  curve 
most  likely  to  be  efficacious  in  an  OP  CWNA  attack.  Most 
oximes  have  a  half-life  of  only  1-2  h  (Voicu  et  al.  2010).  The 
blood  concentration  endpoint  provides  a  measure  of  TMB-4 
being  delivered  to  OP-inhibited  AChE  in  the  tissues.  Reac¬ 
tivation  of  brain  AChE,  specifically  in  the  ponto-medullar 
area,  is  paramount  to  survival  following  an  OP  exposure 
(Voicu  et  al.  2010),  making  the  brain  concentration  of  TMB-4 
the  other  critical  endpoint  of  choice. 

All  physiological,  physicochemical  and  route-specific 
parameters  for  both  intravenous  and  intramuscular  routes 
were  varied  in  the  sensitivity  analysis.  Only  two  parameters 
were  shown  to  result  in  model  sensitivity  across  endpoints 
for  both  routes  of  administration:  volume  of  venous  blood 
(VVC)  and  volume  of  the  slowly  perfused  tissue  com¬ 
partment  (VSC).  Both  parameters  are  taken  from  published 
physiological  values  (Table  2).  The  model  is  sensitive  to 
VVC  only  during  the  first  few  minutes  of  either  exposure 
route.  The  model  assumes  instant  mixing  once  the  TMB-4 
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enters  the  blood  from  either  the  intravenous  syringe  or  the 
IM  compartment;  small  changes  in  the  VVC  result  in  rel¬ 
atively  large  changes  in  concentration  in  the  blood  and 
therefore  in  the  brain.  VSC  represents  the  largest  collective 
tissue  volume  of  the  body  (62.7  %  of  body  weight, 
Table  4),  which  greatly  affects  distribution  of  the  water- 
soluble  TMB-4.  Although  neither  of  these  physiological 
values  was  altered  in  fitting  the  TMB-4  model,  the  tissue/ 
blood  partition  coefficient  for  the  slowly  perfused  tissue 
compartment,  PS,  was  found  to  be  critical  in  fitting  the  rat 
and  human  parameterizations  of  the  model.  The  sensitivity 
of  the  model  output  to  the  slowly  perfused  tissue  com¬ 
partment  in  general  indicates  that  the  sheer  size  of  this 
compartment  and  the  TMB-4  that  diffuses  into  it  deter¬ 
mines  what’s  left  over  for  the  rest  of  the  tissues  in  the 
model,  blood  and  brain  included. 

An  analysis  of  sensitivity  is  not  complete  without  a  look 
at  the  uncertainty  of  the  parameters  to  which  the  endpoints 
are  sensitive.  The  model  was  sensitive  to  several  physio¬ 
logical  parameters;  specifically,  the  model  was  sensitive  to 
VVC  and  VSC  across  administration  routes  and  endpoints. 
Yet,  these  parameters  are  based  on  physiological  mea¬ 
surement,  reducing  their  uncertainty.  Other  physiological 
parameters  would  also  be  considered  to  have  a  low  degree 
of  uncertainty;  however,  physicochemical  parameters  used 
in  this  model,  such  as  maximum  rate  constants  and  tissue/ 
blood  PCs,  are  more  uncertain  because  they  are  either  fit  or 
predicted  utilizing  QSPR  algorithms.  While  all  PCs  used  in 
the  model  would  be  considered  to  have  a  relatively  higher 
degree  of  uncertainty,  the  model  was  found  to  be  sensitive 
to  only  the  brain/blood  PC  (PBr),  the  IM  compartment  tis¬ 
sue/blood  PC  (PIM)  and  the  kidney/blood  PC  (PK).  Other 
uncertain  physicochemical  parameters  to  which  the  model 
is  sensitive  include  the  blood/brain  diffusion  limitation 
coefficient  (PABrC)  and  the  excretion  parameters  (maxi¬ 
mum  saturable  excretion  rate  constant  (VMaxUrC)  and 
excretion  affinity  constant  (KmUr)).  Additional  laboratory 
studies  designed  to  measure  even  a  few  of  these  parameters 
would  reduce  the  overall  uncertainty  in  the  model. 

Model  application  and  future  work 

The  purpose  of  this  TMB-4  model  is  to  act  as  a  stepping 
stone  for  prediction  of  MMB-4  kinetics.  Unfortunately, 
published  MMB-4  pharmacokinetic  data  were  not  found 
for  any  of  the  three  species  for  which  this  model  was 
developed.  However,  with  consideration  for  the  assump¬ 
tions  made,  the  model  can  predict  anticipated  MMB-4 
kinetics  in  humans.  Figure  12  depicts  the  simulated 
Vojvodic  (1970)  TMB-4  study  kinetics  from  Fig.  6a  along 
with  predicted  kinetics  of  an  equimolar  MMB-4  dose  given 
under  identical  conditions.  Aside  from  the  lower  molecular 
weight  of  MMB-4  compared  to  TMB-4  (258.27  vs. 
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286.33  g/mol),  only  predicted  tissue/blood  PC  values 
changed  between  model  parameterizations  (Table  6).  Fit 
tissue/blood  PC  values  (PBr  and  PIM)  and  tissue  perme¬ 
ability  area  cross  products  remained  the  same  as  in  the 
TMB-4  parameterization.  The  model  predicts  a  slightly 
lower  concentration  of  MMB-4  as  compared  to  the  equi¬ 
molar  dose  of  TMB-4,  due  to  the  lesser  molecular  weight 
of  MMB-4.  Since  the  predicted  pharmacokinetics  of  these 
oximes  are  very  similar,  the  improved  equimolar  efficacy 
of  MMB-4  over  TMB-4  (Shih  et  al.  2009)  is  likely  a 
pharmacodynamic  difference. 

A  key  strength  of  PBPK  models  is  that  they  can  be 
utilized  to  predict  human  blood  levels  of  a  drug  over  a 
variety  of  dosing  regimens.  Following  the  comparison 
made  in  Fig.  12,  the  model  was  utilized  to  predict  a  pre¬ 
sumably  efficacious  dose  in  humans  under  a  dosing  regi¬ 
men  consistent  with  military  deployment  (Fig.  13).  Under 
the  current  protocol,  soldiers  are  issued  three  MARK  I  kits 
(Newmark  2009)  and  are  instructed  to  self-  or  buddy- 
administer  one  kit  immediately  at  the  first  signs  of  nerve 
agent  poisoning.  Up  to  2  more  kits  may  be  administered 
every  10-15  min  if  symptoms  persist  (USAF  1996).  Proper 
use  of  autoinjectors  involves  holding  the  unit  in  place  for 
10  s  to  ensure  complete  transfer  of  the  drug  to  the  patient 
(CHEMM  2011).  Using  these  administration  guidelines, 
the  model  was  iteratively  run  using  progressive  doses  until 
a  presumably  efficacious  blood  concentration  lasting  more 
than  2  hours  was  achieved.  In  vitro  data  from  Worek  et  al. 
(2010)  indicate  that  human  erythrocytes  with  AChE  levels 
inhibited  in  excess  of  95  %  by  sarin,  cyclosarin,  VX  or  VR 
may  be  reactivated  to  40  %  AChE  activity  by  a  media 
concentration  of  50  pM  MMB-4,  theoretically  equivalent 
to  a  blood  concentration  of  12.9  mg  MMB-4  cation/L.  The 
model  predicts  that  3  autoinjectors,  each  containing  the 
equivalent  of  5  mg  MMB-4/kg,  would  result  in  a  poten¬ 
tially  efficacious  dose  lasting  more  than  2  h  in  a  70  kg 
person.  This  scenario  assumes  that  the  soldier  would 
receive  medical  assistance  beyond  buddy  care  before  the 
3-h  mark;  additional  autoinjectors  would  be  needed  to 
maintain  an  efficacious  dose  beyond  that  time  point  if 
medical  care  did  not  arrive.  A  significant  concern  with  this 
scenario  is  the  length  of  time  necessary  for  the  blood 
concentration  to  reach  the  efficacy  level  (21  min),  pro¬ 
viding  the  soldier  received  a  CWNA  dose  high  enough  to 
result  in  more  than  95  %  inhibition  of  AChE.  Adminis¬ 
tering  all  three  autoinjectors  one  after  the  other  brings  the 
time  to  efficacy  down  to  9.24  min,  assuming  autoinjectors 
were  given  40  s  apart  to  allow  time  for  proper  package 
handling  (simulation  not  shown).  These  model  predictions 
do  raise  concern  that  the  current  method  of  administering 
oximes  in  the  held  may  not  result  in  an  efficacious  blood 
concentration  rapidly  enough  to  save  the  lives  of  those 
exposed  to  high  levels  of  CWNAs. 
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Fig.  12  Comparison  of  TMB-4  model  simulation  of  Vojvodic  (1970) 
study  protocol  (solid  line,  see  Fig.  6)  and  prediction  of  an  equimolar 
MMB-4  dose  (dotted  line).  Model  conditions  were  set  for  a  81.9  kg 
male  human  receiving  either  2.45  mg  TMB-4/kg  BW  or  2.21  mg 
MMB-4/kg  BW,  dosed  intramuscularly 


Fig.  13  Prediction  (line)  of  presumably  efficacious  dose  of  MMB-4 
for  military  use.  Model  conditions  were  set  for  a  70  kg  human 
administered  3  autoinjectors  10  min  apart  containing  the  equivalent 
of  5  mg  MMB-4/kg.  Open  diamonds  indicate  successive  doses  while 
the  dotted  lines  indicate  the  period  of  time  over  which  the  blood 
concentration  exceeded  the  potential  efficacious  level  of  12.9  mg/L 
per  in  vitro  data  by  Worek  et  al.  (2010) 


Table  6  Tissue/blood  partition  coefficents  used  in  the  MMB-4  and 
TMB-4  model  parameterizations  for  humans 


Tissue/blood  partition  coefficients  (unitlessf 

MMB-4 

TMB-4b 

Brain 

PBr 

0.13* 

0.13# 

Diaphragm 

PD 

0.94 

0.94 

Fat 

PF 

0.24 

0.19 

IM  compartment 

PIM 

2.5* 

2.5# 

Kidney 

PK 

0.94 

0.94 

Liver 

PL 

0.90 

0.90 

Rapidly  perfused 

PR 

0.98 

0.97 

Slowly  perfused 

PS 

0.74 

0.74 

a  PC  values  predicted  unless  otherwise  noted;  b  Values  from  Table  3 
*  Value  from  TMB-4  model;  #Fit  value 


MMB-4  in  human  volunteers  would  help  support  species 
extrapolation  assumptions. 

While  the  MMB-4  parameterization  of  the  model  is  still 
awaiting  further  development  should  the  opportunity  arise 
in  the  form  of  newly  available  kinetic  data,  the  TMB-4 
PBPK  model  parameterization  can  now  be  expanded  to 
incorporate  pharmacodynamic  effects,  including  oxime 
reactivation  of  OP-AChE  complexes,  to  better  extrapolate 
efficacy  across  the  human  population  under  different  dos¬ 
ing  regimens.  The  resulting  PBPK-pharmacodynamic 
mathematical  model  would  allow  in  silico  testing  of 
countermeasure  dosing  regimens  and  facilitate  extrapola¬ 
tion  of  novel  oxime  effects  in  experimental  animals  to 
humans,  thereby  optimizing  therapeutic  efficacy. 


Although  the  model  can  be  used  in  its  present  state  to 
help  predict  MMB-4  kinetics,  further  development  and 
model  validation  is  recommended  utilizing  MMB-4  kinetic 
data.  Since  rats  are  considered  poor  animal  models  of 
CWNA  effects  due  to  the  abundance  of  carboxylesterase  in 
rat  blood  and  its  bioscavenging  activity  against  OP  com¬ 
pounds  (Maxwell  and  Brecht  2001),  further  model  devel¬ 
opment  for  oximes  as  OP  reactivators  should  be  pursued  in 
another  species,  such  as  the  guinea  pig.  Therefore,  guinea 
pig  blood,  liver,  kidney,  brain  and  urine  time  course 
kinetics  of  MMB-4  would  provide  a  solid  basis  for  further 
model  development.  Additionally,  blood,  urine  and  tissue 
kinetics  in  non-human  primates  would  allow  some  insight 
into  the  magnitude  of  species  differences  in  tissue/blood 
PCs,  metabolism  and  urinary  excretion  highlighted  in  the 
current  phase  of  the  model.  Blood  and  urine  kinetics  for 
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